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Abstract 



The production of jets is studied in deep-inelastic e + p scattering at low negative four mo- 
mentum transfer squared 5 < Q 2 < 100 GeV 2 and at inelasticity 0.2 < y < 0.7 using 
data recorded by the HI detector at HERA in the years 1999 and 2000, corresponding to an 
integrated luminosity of 43.5 pb _1 . Inclusive jet, 2-jet and 3-jet cross sections as well as 
the ratio of 3-jet to 2-jet cross sections are measured as a function of Q 2 and jet transverse 
momentum. The 2-jet cross section is also measured as a function of the proton momentum 
fraction £. The measurements are well described by perturbative quantum chromodynam- 
ics at next-to-leading order corrected for hadronisation effects and are subsequently used to 
extract the strong coupling a s . 
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1 Introduction 



Jet production in neutral current (NC) deep-inelastic scattering (DIS) at HERA provides an 
important testing ground for Quantum Chromodynamics (QCD). While inclusive DIS gives 
only indirect information on the strong coupling via scaling violations of the proton structure 
functions, the production of jets allows a direct measurement of a s . The Born level contribution 
to DIS (Fig.QJa)) generates no transverse momentum in the Breit frame, where the virtual boson 
and the proton collide head on [OQ. Significant transverse momentum P T in the Breit frame is 
produced at leading order (LO) in the strong coupling a s by the QCD Compton (Fig. (Ub)) 
and boson-gluon fusion (Fig. HJc)) processes. The latter dominates jet production for the range 
of the negative four momentum transfer squared of this analysis, 5 < Q 2 < 100 GeV 2 , and 
provides direct sensitivity to the gluon density function of the proton [|2]|. 

Analyses of inclusive and multi-jet production in DIS were previously performed at high Q 2 
(> 100 GeV 2 ) AMI and at low Q 2 (< 100 GeV 2 ) MM by the HI and ZEUS Collaborations 
at HERA. In this paper new measurements of the inclusive jet, 2-jet and 3-jet production cross 
sections, as well as the ratio of 3-jet to 2-jet cross sections, are presented as a function of Q 2 and 
the jet transverse momenta in the Breit frame, P T , in the ranges 5 < Q 2 < 100 GeV 2 and P T > 
5 GeV. The 2-jet cross section is also presented as a function of £ = £Bj(l + M 2 2 /Q 2 ), which 
in LO corresponds to the momentum fraction of the proton carried by the interacting parton (see 
Figsffib) andQJc)). The variable xbj denotes the Bjorken scaling variable and Mi 2 the invariant 
mass of the two jets of highest Pt- The data correspond to higher integrated luminosity and a 
higher centre-of-mass energy than in the previous HI analyses at low Q 2 (2l[6l|. The larger data 
set together with improved understanding of the hadronic energy measurement significantly 
reduces the total uncertainty of the cross section measurements. The results are compared 
with perturbative QCD predictions at next-to-leading order (NLO) corrected for hadronisation 
effects, and a s is extracted from a fit of the predictions to the data. These measurements allow 
the running of the strong coupling to be tested down to the limits of the perturbative calculation. 
Together with the high Q 2 measurements flU the data test the running of a s in the range of 
renormalisation scale ji r between about 6 and 70 GeV. 



Figure 1: Deep-inelastic lepton-proton scattering at different orders in a s : (a) Born contribution 
0(1) and 0(a s ) processes (b) QCD Compton scattering and (c) boson-gluon fusion. 
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2 Experimental Method 



The data used for this analysis were recorded with HI detector in the years 1999 and 2000, when 
HERA collided positrons of energy E e = 27.6 GeV with protons of energy E p = 920 GeV 
giving a centre-of-mass energy \/s = 319 GeV. The corresponding integrated luminosity is 
43.5 pb" 1 . 

2.1 The HI detector 

A detailed description of the HI detector can be found in EUl. Here, a brief account of the 
components most relevant to the present analysis is given. The origin of the HI coordinate 
system is the nominal ep interaction point. The direction of the proton beam defines the positive 
2-axis (forward direction). The polar angle 9 is measured with respect to this direction. The 
pseudorapidity is defined as rj = — lntan(0/2). 

In the central region (20°<#<160°) the ep interaction region is surrounded by a two- 
layered silicon strip detector [fT0| and two large concentric drift chambers (CJCs), operated in- 
side a 1.16 T solenoidal magnetic field. The trajectories of charged particles are measured in the 
central tracker with a transverse momentum resolution of ct(pt)/pt = 0.006 px/GeV © 0.02 |fi~Tfl 
Two additional drift chambers complement the CJCs by precisely measuring the ^-coordinates 
of track segments and hence improve the determination of the polar angle. The central tracking 
detectors also provide triggering information based on track segments measured in the r-<f> plane 
of the central jet chambers and on the z position of the event vertex obtained from the double 
layers of two multiwire proportional chambers (MWPCs). The forward tracking detector and 
the backward drift chamber (BDC) measure tracks of charged particles at smaller (7° <9< 25°) 
and larger (155° < 9 < 175°) polar angle than the central tracker, respectively. 

A finely segmented electromagnetic and hadronic liquid argon (LAr) calorimeter lfi~2l sur- 
rounds the tracking chambers. It has a polar angle coverage of 4°<#<154° and full az- 
imuthal acceptance. The energy resolution is a (E)/E = 0.12/^/ GeV © 0.01 for elec- 
tromagnetic showers and a{E)/E = 0.5/ \JE~/ GeV © 0.02 for hadrons, as measured in test 
beams lfi~3l . A lead-scintillating fibre spaghetti calorimeter (SpaCal) flU covers the backward 
region 153°<#<178°. Its main purpose is the detection of scattered positrons. The energy 
resolution of the SpaCal for positrons is a(E) /E = 0.071/ y/E/ GeV © 0.01. 

The luminosity is measured via the Bethe-Heitler Bremsstrahlung process ep — > ep"f, the 
final state photon being detected in a crystal calorimeter at z = — 103 m. 

2.2 Event and jet selection 

The data sample of this analysis was collected using a combination of triggers which require the 
scattered positron to be measured in the SpaCal, at least one high transverse momentum track 
(pr > 800 MeV) to be reconstructed in the central tracking chambers and an event vertex to be 
identified by the MWPCs. The trigger efficiency is close to 100% for the whole analysis phase 
space as determined from the data using independent triggers as a reference. 
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The selection of NC DIS events is based on the identification of the scattered positron as the 
most energetic compact calorimetric deposit in the SpaCal with an energy E' e > 7.5 GeV and 
a polar angle 156° < 9' e < 175°. The energy weighted radius of this cluster is required to be 
less than 4 cm, as expected for an electromagnetic shower. The cluster must be geometrically 
associated with a track candidate in the BDC. The ^-coordinate of the primary event vertex is 
required to be within ±35 cm of the nominal position of the interaction point. 

The remaining clusters in the calorimeters and the charged tracks are combined to re- 
construct the hadronic final state, using an algorithm which avoids double counting of en- 
ergy |fi~4l[T5ll . The total longitudinal energy balance, determined as the difference of the total 
energy E and the longitudinal component of the total momentum P z , calculated from all de- 
tected particles including the scattered positron, must satisfy 45 < E — P z < 65 GeV. This 
requirement reduces contributions of DIS events with hard initial state photon radiation. For 
the latter events, the undetected photons propagating in the negative z direction lead to values 
of E — P z significantly lower than the expected value 2E e = 55.2 GeV. After this selection the 
contribution from photoproduction is negligible as estimated using Monte Carlo simulations. 

The kinematic region covered by this analysis is defined by 

5<Q 2 <100GeV 2 and 0.2<y<0.7, 

where y = Q 2 /(s ■ xb; } ) quantifies the inelasticity of the interaction. These two variables are 
reconstructed from the four momenta of the scattered positron and the hadronic final state par- 
ticles using the electron-sigma method IfToll . 

Jet finding is performed in the Breit frame. The boost from the laboratory system is deter- 
mined by Q 2 , y and by the azimuthal angle of the scattered positron. Particles of the hadronic 
final state are clustered into jets using the inclusive kx algorithm lfi~7l[T8l with the massless Pt 
recombination scheme and with the distance parameter in the r] — <p plane R = 1. The cut 
— 1 < r/Lab < 2-5, where r]^ ah is the jet pseudorapidity in the laboratory frame, ensures that jets 
are contained within the acceptance of the LAr calorimeter. The transverse energy of jets in the 
Breit frame is required to be above 5 GeV. Jets are ordered by decreasing transverse momentum 
P T in the Breit frame, which is identical to the transverse energy E T for massless jets. The jet 
with highest Pt is referred to as the "leading jet". 

Three jet samples are defined: the inclusive jet sample contains all jets which satisfy the jet 
selection criteria; the 2-jet and 3-jet samples contain events with at least 2 and 3 jets, respec- 
tively. In addition, to avoid regions of phase space where fixed order perturbation theory is not 
reliable |fi~9l , 2-jet events are accepted only if the invariant mass M\i of the two leading jets 
exceeds 18 GeV. The same requirement is applied for 3-jet events such that the 3-jet sample is 
a subset of the 2-jet sample. 

The selection criteria are summarised in Table [Q The final inclusive jet sample contains 
164522 events with 230140 jets. The 2-jet and 3-jet samples contain 31550 and 4879 events 
respectively. 
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2.3 Determination of the jet cross sections 

In order to extract the jet cross sections at hadron level, the experimental data are corrected bin- 
by-bin for effects of limited detector acceptance and resolution and for QED radiation effects. 
The following leading order Monte Carlo event generators are used for the correction procedure: 
DJANGOH EOl . which uses the Colour Dipole Model with QCD matrix element corrections as 
implemented in ARIADNE OTTl . and RAPGAP (23, based on QCD matrix elements matched 
with parton showers in leading log approximation. The effects of QED radiation are included 
using the HERACLES flU program interfaced with RAPGAP and DJANGOH. In both Monte 
Carlo generators the hadronisation is modelled with Lund string fragmentation [|24]|. The gen- 
erated events are passed through a GEANT3 [25] based simulation of the HI apparatus and are 
reconstructed using the same program chain as for the data. Both the RAPGAP and DJANGOH 
simulations provide a good overall description of the shapes of all relevant data distributions. 
To further improve the agreement between Monte Carlo and data, a reweighting as a function 
of Q 2 and P T of the leading jet is applied to the Monte Carlo events. 

The bin dependent correction factors are determined from Monte Carlo simulations as the 
ratios of the cross sections obtained from particles at hadron level without QED radiation to the 
cross section calculated using reconstructed particles and including QED radiation effects. The 
mean values of the correction factors determined by RAPGAP and DJANGOH are used, and 
half of the difference is assigned as a model uncertainty. The typical value of these factors is 
between 1.2 and 1.4. 

The binnings in Q 2 , Pt and £ used to measure the jet cross sections are given in Table[2] The 
bin purities, defined as the fraction of events reconstructed in a particular bin that originate from 
that bin on hadron level, is found to be typically 70% and larger than 50% in all analysis bins. 
The bin stabilities, defined as the fraction of events originated from a particular bin on hadron 
level that are reconstructed in that bin, is typically 60% and larger than 40% in all analysis bins. 

2.4 Experimental uncertainties 

Several sources of experimental uncertainties are considered. The systematic uncertainties of 
the jet cross sections are determined by propagating the corresponding estimated measurement 
errors through the full analysis: 

• The relative uncertainty of the positron energy calibration is better than 1%. The absolute 
uncertainty of the positron polar angle is about 1 mrad. Uncertainties in the positron 
reconstruction affect the event kinematics and thus the boost to the Breit frame. This 
in turn leads to a relative error of up to 2% on the jet cross section for each of the two 
sources. 

• The relative uncertainty on the energy of the reconstructed hadronic final state as well as 
of the jet energy is estimated to be 2%. It is dominated by the uncertainty of the hadronic 
energy scale of the LAr calorimeter. Two different calibration methods are used for jet 
transverse momentum below and above 10 GeV, respectively. The resulting uncertainty 
on the cross sections is typically in the range of 4% to 10%. The uncertainty of the 
SpaCal hadronic energy scale of 7% contributes less than 1% to the uncertainty of the 
cross section. 
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• The model dependence of the detector correction factors is estimated as described in 
section [23] It reflects the sensitivity of the detector simulation to the details of the model, 
especially the parton showering and its impact on the migration between adjacent bins in 
Pp. The model dependence is below 10% in most of the bins and typically 4%. 

• The luminosity measurement uncertainty leads to an overall normalisation error of the jet 
cross sections of 1.5%. 

The uncertainty of the luminosity measurement is assumed to be fully correlated between 
the bins. The remaining sources of systematics, namely the positron energy scale and polar 
angle, the hadronic final state energy scale and the model dependence are assumed to be equally 
shared between correlated and uncorrected parts. 

The dominant experimental uncertainties on the jet cross sections arise from the model 
dependence of the data correction and from the LAr hadronic energy scale uncertainty. The 
individual contributions are added in quadrature to obtain the total systematic uncertainty. 



3 NLO QCD Calculations 

Reliable quantitative predictions of jet cross sections in DIS require the perturbative calculations 
to be performed at least to NLO in the strong coupling. The NLO calculations are used for 
comparison to data and for the a s extraction. By using the inclusive k T jet algorithm, the 
observables in the present analysis are infrared and collinear safe. Application of this algorithm 
in the Breit frame allows the initial state singularities to be absorbed in the definition of the 
proton parton densities, as needed for the calculation of factorised jet cross sections ll26l . 

Jet cross sections are predicted at the parton level using the NLOJET++ program 11271 at 
NLO in the strong coupling using the same jet definition as in the data analysis. When com- 
paring data and theory predictions the strong coupling is taken to be a s (M z ) = 0.118 at the 
Z° boson mass and is evolved as a function of the renormalisation scale with two loop preci- 
sion. No QED radiation is included in the calculations, but the running of the electromagnetic 
coupling with Q 2 is taken into account in the theoretical predictions. The calculations are per- 
formed in the MS scheme GBl for five massless quark flavours. The parton density functions 
(PDFs) of the proton are taken from the CTEQ6.5M set ||29ll . The factorisation scale /!/ and 
the renormalisation scale \i r are taken to be a/ (Q 2 + Pp bs)/2 f° r me NLO predictions, with 
-fr.obs denoting the P T of the jet for inclusive jet cross sections and the average transverse mo- 
mentum of the two leading jets (Pt) for the 2-jet and 3-jet cross sections. This choice of the 
scales is motivated by the presence of two hard scales in jet production in DIS, Pt and Q, the 
latter being smaller in most of the analysis bins. The calculations were also performed using 
H r = P;r,obs- With this choice of renormalisation scale the NLO QCD prediction decreases by 
10 — 20% at lowest Q 2 and P T and is disfavoured by the data. 

Hadronisation corrections are calculated for each bin using Monte Carlo event generators 
DJANGOH and RAPGAP which implement different models for parton showering. These cor- 
rections are determined as the ratio of the cross section at hadron level to the cross section at the 
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partem level after parton showers. It was verified that the parton level jet cross sections obtained 
from DJANGOH and RAPGAP are in agreement with those from the NLO calculation within 
the systematic uncertainties considered here. The hadronisation correction factors are deter- 
mined as the average of values obtained from DJANGOH and RAPGAP. Half of the difference 
is assigned as hadronisation uncertainty and included as a part of the theoretical uncertainty. 
For inclusive and 2-jet cross sections the hadronisation correction factors differ typically by less 
than 10% from unity and agree at the level of 2 to 5% between the two Monte Carlo simulations. 
For 3-jet cross sections, as well as for the ratios of 3-jet to 2-jet cross sections, the hadronisation 
correction factors differ from unity by about 20% with up to 10% difference between the two 
MC models. 

The dominant theoretical uncertainty is related to the missing higher orders in the perturba- 
tive calculation, and is conventionally estimated by separately varying the scales /i / and \i T by 
factors in the arbitrary range 0.5 to 2. The contributions from the two scale variations are similar 
and are added in quadrature to obtain the total scale dependence uncertainty. The uncertainty 
originating from the PDFs is taken into account for the a s extraction using the variations of the 
CTEQ6.5M set of parton densities. 



4 Cross Section Measurements 

In the following, the differential cross sections, corresponding to the phase space given in Ta- 
ble [Q are presented for inclusive jet, 2-jet and 3-jet production at hadron level. Ratios of 3-jet 
to 2-jet hadron level cross sections are also presented. The measurements are shown in Tables 
[3] to|§] and Figs. [2] to [8] 

4.1 Inclusive jet cross section 

The measured inclusive jet cross section, corrected for detector and radiative QED effects, is 
presented as a function of Q 2 and Pt of the jet, as single differential distributions in Figs.[2a,b) 
and double differentially in Fig. [3] Each jet which satisfies the jet selection criteria described in 
section [2T21 enters these distributions. 

The measurements are well described by NLO QCD predictions corrected for hadronisation 
effects as explained in section |3] The theoretical uncertainty, dominated by the scale variation, 
reaches 30% for the lowest Q 2 and Pt bins and decreases to 10% for the highest Q 2 and Pt 
values. The relative contribution of hadronisation corrections to this error is small. The PDF 
uncertainty is about 6% at the lowest Q 2 and Pt and decreases to 2% for the highest Q 2 . 

4.2 2-jet and 3-jet cross sections 

The measured single differential cross sections for 2-jet and 3-jet production as functions of Q 2 
and of the average transverse momentum of the two leading jets (Pt) are shown in Figs.[2£c,e) 
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and Figs.[2d,f) and are well described by the NLO QCD calculations corrected for hadronisa- 
tion. The relative uncertainties on the NLO QCD calculations of 2-jet, 3-jet and inclusive cross 
sections are all of similar size. 

The double differential 2-jet cross sections are presented in seven Q 2 bins as functions of 
the variables (Pt) and £ in Fig. |4] and Fig. [5] respectively. The 3-jet cross section is shown in 
four Q 2 bins as a function of (P T ) in Fig. [6j The NLO QCD calculation provides an overall 
good description of the measured distributions within the quoted theoretical and experimental 
uncertainties. Requirements in Pt and M\i suppress the cross section at low £ where a rise is 
expected due to the increase of the gluon density. 

The present results for the inclusive and the 2-jet cross sections were compared to the pre- 
vious HI results in [|2]| and J6]|, respectively. Taking into account the difference between proton 
beam energies and differences in the kinematic region studied, the results are found to be con- 
sistent with each other. 

The 3-jet cross section normalised to the 2-jet cross section is presented in Fig.[7]for single 
differential and in Fig. [8] for double differential distributions. This observable benefits from 
cancellation of the normalisation uncertainties and reduction of the other systematic uncertain- 
ties by about 50%. It is described by the NLO cross section except for the lowest {Pt) bin, 
as seen in Fig. [7£b) and Fig. [81 and shows a reduced sensitivity to the renormalisation scale 
variation which is done simultaneously for 2-jet and 3-jet cross sections. The ratio shows no 
significant dependence on Q 2 (Fig. |7ja)), but increases with (Pt) (Figs. [7£b)|8]) mainly due to 
the increasing phase space. 



5 Extraction of the Strong Coupling 

The QCD predictions for jet production depend on a s and on the PDFs of the proton. The strong 
coupling a s is determined from the measured jet cross sections using the PDFs as obtained from 
inclusive DIS data and other measurements. 

5.1 Data and QCD predictions 

The a s determination is performed from individual bins of the double differential inclusive jet 
cross section, d 2 a/dQ 2 dP T , and the 2-jet and 3-jet cross sections, d 2 a jdQ 2 d (Pt)- Only bins 
are used in which the size of the A;-factor, defined as the ratio of the cross sections calculated 
in NLO and LO (both obtained with NLOJET++), is below 2.5. The other bins most likely are 
affected by slow convergence of perturbation series and exhibit a high scale dependence, up to 
30% at NLO. 

The requirement that the fc-factor be less than 2.5 corresponds to removing all points with 
-fr.obs < 10 GeV for Q 2 < 20 GeV 2 from the inclusive jet and the 2-jet cross sections and 
points with 10 < Pt < 15 GeV for Q 2 < 10 GeV 2 from the inclusive jet cross section Q. In 

'The bins removed from the a s analysis also correspond to the energy regime which is close to the 6-quark 
mass threshold, where the five flavour massless approximation used in NLOJET++ is not expected to be valid. 
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total, 62 cross sections measurements are used for a s extraction: 22 inclusive jet, 24 2-jet and 
all 16 3-jet points. 

QCD predictions of the jet cross sections are calculated as a function of a s (fj, r ) with the 
FastNLO package ll30l using the CTEQ6.5M proton PDFs and applying the hadronisation cor- 
rections as described in section [3] 

5.2 The x 2 definition 

Measurements and theory predictions are used to calculate a x 2 ( a s) with the Hessian method OTTl . 
where parameters representing the systematic shifts of detector related observables, described 
in section [2741 are left free in the fit. The shifts found by the fit are consistent with the a priori 
estimated experimental uncertainties. Due to different calibration strategies for jets with P T 
above and below 10 GeV, two different parameters are used for the hadronic final state energy 
scale for bins with Pr,obs < 10 GeV and Pr, bs > 10 GeV. The Hessian method used here takes 
into account correlations of experimental uncertainties and has also been used in global data 
analyses D 111321 and in previous HI publications 0113311 . The statistical correlations among the 
different bins and different observables are treated as described in flU[34]|. The experimental 
uncertainty of a s is defined by the change in a s which gives an increase in x 2 of one unit with 
respect to the minimal value. 

5.3 Theory and PDF uncertainties 

The theoretical uncertainty is estimated by the offset method as the difference between the value 
of a s from the nominal fit to the value when the fit is repeated with independent variations of 
different sources of theoretical uncertainties as described in section [3] The resulting uncertain- 
ties due to the different sources are summed in quadrature. The up (or down) variations are 
applied simultaneously to all bins in the fit. The impact of hadronisation corrections and the 
factorisation scale uncertainty on a s typically amounts to 1% to 2% for each source. The largest 
uncertainty, of typically 8%, corresponds to the accuracy of the NLO approximation to the jet 
cross section estimated by varying the renormalisation scale as described in section [3l 

The uncertainty due to PDFs is estimated by propagating the CTEQ6.5M errors. The typical 
size of the resulting error is 2% for a s determined from the inclusive jet or 2-jet cross sections 
and 1% when measured with the 3-jet cross sections. This uncertainty is twice as large as that 
estimated with the uncertainties given for the MSTW2008nlo90cl set [|35l which in turn exceeds 
the difference between a s values extracted with the central sets of CTEQ6.5M, CTEQ6.6M ll36l 
and MSTW2008nlo. 

The CTEQ6.5M PDF parameterisation was obtained assuming a s (M z ) = 0.118. In order 
to test whether this value of a s (Mz) biases the results obtained using the nominal method 
presented above, a method, similar to the one used in 01, is employed using the PDFs from the 
CTEQ6.6 series, which were obtained assuming different values for a s (Mz). The cross section 
as a function of the strong coupling is interpolated with a polynomial and this interpolation 
is used to determine the best fit of the strong coupling to the data. The result obtained with 
this alternative fit method is found to be compatible, well inside one standard deviation of 
the experimental error, with the value determined by the nominal method. Hence there is no 
indication for a bias due to the value of the strong coupling assumed for the CTEQ6.5M PDFs. 
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5.4 Fit results 



The fits of the strong coupling a s are performed individually for each of the 62 cross section 
measurements as described in section l5Tl These measurements constrain the value of the strong 
coupling at the Z° mass, a s (Mz). As an example, ot s {Mz) values determined from the inclusive 
jet cross section are shown in Fig. [9l Different a s (Mz) values agree within experimental errors. 
For each of the three observables and each Q 2 region, fits of a s (M z ) to all Pt or (P T ) bins 
in that region are performed. The fit results are evolved from Mz to the average scale jj r in 
the respective Q 2 regions, where the average renormalisation scale fi r is calculated using NLO 
predictions. Figs. [TOta-c) show the obtained a s (fi r ) values for fits to inclusive, 2-jet and 3-jet 
cross-sections, respectively. Also shown in these figures are QCD predictions a s (fi r ), derived 
from common fits of a s (M z ) to all respective Q 2 and Pt or (P T ) bins. The results of these 
three common fits are summarised in Table [TOl If the points with fc-factor above 2.5 are also 
included, the a s (Mz) values obtained are changed by less than one standard deviation of the 
total experimental uncertainty. As the a s (Mz) measurements derived from inclusive, 2-jet or 
3-jet observables agree within uncertainties in any of the Q 2 regions, they are combined within 
four Q 2 regions, taking into account statistical and experimental systematic correlations. These 
results, evolved from the scale M z to the average \i r in each region, are shown in Fig. [TIT a). 
As compared to the case of extracting a s from only one variable, the experimental uncertainties 
are reduced significantly. 

Finally, all 62 data points are used in a common fit of the strong coupling taking the corre- 
lations into account with a fit quality x 2 /ndf = 49.8/61: 

a s (M z ) = 0.1160 ± O.OOU{exp.)toSi th -) ± 0.0016(pdf). (1) 

The experimental error on a s (M z ) measured with each observable typically amounts to 1.5%. 
The combination of different observables, even though partially correlated, gives rise to ad- 
ditional constraints on the strong coupling and leads to an improved experimental uncertainty 
of 1.2%. This error changes by at most 20% when the scales, hadronisation factors and PDF 
parameterisation are changed within the limits defined by their uncertainties. The total error is 
dominated by the theoretical uncertainty of about 7% mainly due to scale variations. 

The determination of the strong coupling from the ratio of the 3-jet to the 2-jet cross section 
provides an alternative approach to combining the different cross section data. On the one hand 
the sensitivity of this observable to a s , which is 0(a s ), is reduced with respect to the 3-jet cross 
section, which is 0(a 2 ). On the other hand this observable benefits from reduced experimental 
and theoretical uncertainties (see section |4T2l) . The common fit of the strong coupling to the 14 
ratio points, for which the ^-factors are below 2.5 for both the 3-jet and 2-jet cross sections, is 
given in Table [TOl The extracted a s as a function of fi r is shown in Fig.QjJb). The experimental 
uncertainty on a s increases to 3% with respect to the combined fit from the cross sections. The 
theoretical uncertainty is reduced to 5% and is dominated by the hadronisation uncertainty. 

The strong coupling extracted from all 62 data points (Q~|) agrees well with that obtained from 
jet cross sections in the higher Q 2 range between 150 and 15000 GeV 2 : a s (M z ) = 0.1168 ± 
0.0007(exp.)iooo3o(^-) ±0.0016(pdf) 0). This agreement is remarkable, given the sensitivity 
of the NLO prediction to the renormalisation scale in the low Q 2 regime. At high Q 2 flU 
the fit was done using the factorisation scale Q, instead of a/ (Q 2 + P| j0bs )/2 used in this 
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analysis. However, as was shown in flU, the choice of factorisation scale has only little impact on 
extracted value of a s at high Q 2 . The value of a s (M z ) obtained in this analysis is also consistent 
with the world averages a s {M z ) = 0.1176±0.0020 (33 mda s (M z ) = 0.1184 ±0.0007 11581. 

The new low Q 2 measurement together with data from the high Q 2 analysis provides a test 
of the running of the strong coupling for fi r between 6 and 70 GeV as illustrated in Fig. [121 
A simultaneous fit of a s from low and high Q 2 data was also performed. It did not lead to an 
improved precision with respect to the high Q 2 determination alone due to the large theoretical 
uncertainties at low Q 2 . 



6 Conclusion 

Measurements of the inclusive, 2-jet and 3-jet cross sections in deep-inelastic positron-proton 
scattering are presented in the range 5 < Q 2 < 100 GeV 2 and 0.2 < y < 0.7. Jets are recon- 
structed using the inclusive kx algorithm in the Breit frame and are required to have a minimum 
transverse momentum of 5 GeV. Calculations at NLO QCD, corrected for hadronisation effects, 
provide a good description of the single and double differential cross sections as functions of 
the jet transverse momentum P T , the boson virtuality Q 2 as well as of the proton momentum 
fraction £. The precision of the measurements is typically 6 to 10%. 

The strong coupling a s is determined from a fit of the NLO prediction to the measured 
jet cross sections. The dominant source of uncertainties is related to the renormalisation scale 
dependence, which is used to estimate the effect of missing higher orders. The extracted value 
of the strong coupling 

a s (M z ) = 0.1160 ± 0.0014(ezp.)l£^(t/i.) ± 0.0016(pdf) 

is consistent with the value determined from high Q 2 jet cross sections. Both measurements test 
a running of the strong coupling for renormalisation scales \x r between 6 and 70 GeV. 
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NC DIS Selection 


5 < Q 2 < 100 GeV 2 , 0.2 < y < 0.7 


Inclusive jet 


P T > 5 GeV 




-1-0 < <1 < 2.5 


2-jet 


If- 1 , P$ t2 >5GeV 


M 12 > 18 GeV 


3-jet 


Pl ctl , if t2 , pj? t3 > 5 GeV 



Table 1 : The NC DIS and jet selection criteria. 



bin number 
(inclusive and 2-jets) 


corresponding Q 2 range 


1 
2 
3 
4 

J 

6 
7 


5 < Q 2 < 7 GeV 2 
7 < Q 2 < 10 GcV 2 
10 < Q 2 < 15 GcV 2 
15 < Q 2 < 20 GcV 2 

on o2 on ^\T^ 

30 < Q 2 < 40 GcV 2 
40 < Q 2 < 100 GeV 2 


bin number 
(3-jets and 3-jets/2-jets) 


corresponding Q 2 range 


I 

II 

TTT 
111 

IV 


5 < Q 2 < 10 GcV 2 
10 < Q 2 < 20 GcV 2 
on ^ o 2 ^ /in rw 2 

40 < Q 2 < 100 GcV 2 


bin letter 


corresponding Pt or (Pt) range 


a 
b 
c 
d 


5 < P T < 10 GeV 
10 < P T < 15 GeV 
15 < P T < 20 GcV 
20 < P T < 80 GcV 


bin letter 


corresponding £ range 


A 
B 
C 
D 

E 
F 


0.004 < £ < 0.006 
0.006 < £ < 0.010 
0.010 < £ < 0.025 
0.025 < £ < 0.050 
0.05 < £ < 0.1 
0.1 < £ < 0.3 



Table 2: Nomenclature for the bins in negative four momentum transfer squared Q 2 , jet trans- 
verse momentum P T , average transverse momentum of the two leading jets (Pt) and momen- 
tum fraction £ used in the following tables. 
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Inclusive Jet Cross Section ^M- 











total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 






statistics 




uncorrected 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




section 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




[pb/GeV 2 l 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 


633 


0.7 


6.7 


4.7 


4.8 


1.7 


0.3 


0.4 


4.2 


0.88 


5.8 


2 


421 


0.6 


7.0 


4.9 


5.1 


2.2 


0.2 


0.4 


4.3 


0.89 


5.3 


3 


250 


0.6 


6.7 


4.7 


4.9 


1.4 


0.1 


0.7 


4.4 


0.89 


4.7 


4 


158 


0.7 


6.6 


4.6 


4.8 


1.4 


0.1 


0.7 


4.3 


0.90 


4.0 


5 


96.9 


0.7 


6.5 


4.5 


4.7 


1.3 


0.2 


0.5 


4.2 


0.91 


3.5 


6 


59.7 


0.9 


6.4 


4.4 


4.6 


1.5 


0.1 


0.5 


4.0 


0.91 


2.6 


7 


21.0 


0.7 


5.9 


4.1 


4.3 


1.3 


0.1 


0.6 


3.7 


0.92 


1.5 










2-Jet Cross Section ^ 

dQ 


jet 
2 
















total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


cross 


statistical 


total 


uncorrected 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




section 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




[pb/GeV 2 ] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 


83.6 


1.5 


6.8 


4.8 


4.8 


1.2 


0.2 


0.2 


4.4 


0.91 


2.8 


2 


56.4 


1.4 


8.4 


5.9 


6.0 


3.6 


0.1 


0.5 


4.5 


0.93 


2.9 


3 


33.9 


1.3 


7.4 


5.2 


5.2 


1.5 


0.1 


0.7 


4.7 


0.93 


1.8 


4 


22.5 


1.6 


7.7 


5.5 


5.4 


2.0 


0.1 


0.7 


4.8 


0.94 


1.8 


5 


14.4 


1.4 


7.2 


5.1 


5.1 


1.5 


0.2 


0.5 


4.6 


0.94 


2.0 


6 


9.70 


1.9 


7.3 


5.2 


5.1 


1.5 


0.1 


0.6 


4.6 


0.94 


1.9 


7 


3.72 


1.4 


7.2 


5.1 


5.1 


1.2 


0.2 


0.6 


4.7 


0.95 


1.9 



3-Jet Cross Section 











total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


cross 


statistical 


total 


uncorrected 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




section 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




[pb/GeV 2 ] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 


10.5 


2.6 


10.8 


7.7 


7.5 


4.6 


0.4 


0.4 


5.6 


0.81 


6.3 


II 


4.56 


2.7 


10.6 


7.7 


7.4 


4.0 


0.7 


0.8 


5.9 


0.82 


5.5 


III 


2.05 


2.9 


10.8 


7.9 


7.5 


4.1 


0.7 


0.6 


6.0 


0.80 


7.2 


IV 


0.560 


3.8 


11.7 


8.6 


7.9 


4.5 


1.0 


0.6 


6.2 


0.81 


6.1 



Table 3: Inclusive jet, 2-jet and 3-jet cross sections in NC DIS measured as a function of Q~. 
The measurements refer to the phase-space given in table [Q In the columns 3 to 6 are shown 
the statistical uncertainty, the total experimental uncertainty, the total uncorrected uncertainty 
including the statistical one and the total bin-to-bin correlated uncertainty calculated as the 
quadratic sum of the following five components: the model dependence, the electron and the 
hadronic energy scales, the electron polar angle and the luminosity measurement uncertainties. 
The contributions to correlated uncertainties are listed in columns 7 to 10. The sharing of the 
uncertainties between correlated and uncorrelated sources is described in detail in section l2~4l 
The hadronisation correction factors applied to the NLO predictions and their uncertainties are 
shown in columns 1 1 and 12. The bin nomenclature of column 1 is defined in table [2] 
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Inclusive Jet Cross Section 











total 


total 


single contributions to correlated uncertainty 




hadronisation 


hadronisation 


bin 


cross 
section 

[pb/GeV] 


statistical 
uncert. 
[%] 


total 
uncert. 
[%] 


uncorrelated 
uncertainty 
[%] 


correlated 
uncert. 
[%] 


model 
uncert. 
[%] 


electron 
energy scale 
[%] 


electron 
polar angle 
[%] 


hadronic 
energy scale 
[%] 


correction 
factor 


correction 
uncertainty 
[%] 


a 
b 
c 
d 


1250 
174 
38.3 
1.53 


0.3 
0.7 
1.5 
2.2 


6.8 
6.6 

9.8 
15.3 


4.7 
4.5 
6.9 
10.9 


4.9 
4.7 
6.9 
10.8 


2.1 
2.1 
3.4 
8.1 


0.2 
0.1 

0.3 
0.6 


0.4 
0.5 
0.6 
0.8 


4.1 
3.9 
5.8 
6.9 


0.89 
0.93 
0.95 
0.95 


4.4 
2.8 
2.2 
1.5 










2-Jet Cross Section 


do"2-jct 
d(P T > 




















total 


total 


single contributions to correlated uncertainty 




hadronisation 


hadronisation 


bin 


cross 
section 

[pb/GeV] 


statistical 
uncert. 
[%] 


total 
uncert. 
[%] 


uncorrelated 
uncertainty 
[%] 


correlated 
uncert. 
[%] 


model 
uncert. 
[%] 


electron 
energy scale 
[%] 


electron 
polar angle 
[%] 


hadronic 
energy scale 
[%] 


correction 
factor 


correction 
uncertainty 
[%] 


a 
b 
c 
d 


120 
72.5 
15.5 
0.620 


0.8 
0.9 
2.0 
3.0 


6.8 
7.7 
11.1 
17.8 


4.7 
5.4 
7.9 
12.7 


4.9 
5.5 
7.8 
12.5 


2.1 
2.3 
4.9 
9.7 


0.2 
0.1 
0.4 
0.6 


0.4 
0.5 
0.6 
0.8 


4.1 
4.7 
5.9 
7.7 


0.92 
0.94 
0.95 
0.96 


2.1 
2.3 
2.3 
1.8 



3-Jet Cross Section 











total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


cross 


statistical 


total 


uncorrelated 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




section 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




[pb/GeV] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


a 


11.7 


2.5 


9.2 


6.6 


6.3 


2.4 


0.5 


0.5 


5.6 


0.77 


9.5 


b 


16.1 


2.1 


9.0 


6.4 


6.3 


2.5 


0.6 


0.6 


5.5 


0.81 


6.5 


c 


4.29 


3.9 


12.7 


9.4 


8.6 


5.3 


0.8 


0.7 


6.6 


0.84 


3.6 


d 


0.192 


5.6 


20.0 


14.7 


13.6 


10.6 


1.0 


0.9 


8.4 


0.87 


4.1 



Table 4: Inclusive jet, 2-jet and 3-jet cross sections in NC DIS measured as a function of Pt 
for inclusive jet and (Pt) of the two leading jets for 2-jet and 3-jet cross sections together with 
their relative errors and hadronisation correction factors. Other details are given in the caption 
to table The bin nomenclature is defined in tabled 
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Inclusive Jet Cross Section 











total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


cross 


statistical 


total 


uncorrelated 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




section 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




[pb/GoV 3 ] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 a 


109 


0.7 


7.1 


4.9 


5.1 


2.1 


0.4 


0.5 


4.3 


0.87 


6.7 


1 b 


12.9 


1.8 


6.1 


4.4 


4.2 


2.7 


0.1 


0.1 


2.9 


0.91 


3.2 


1 c 


2.77 


3.8 


10.1 


7.6 


6.7 


3.4 


0.2 


0.1 


5.6 


0.95 


3.4 


1 d 


0.104 


6.4 


15.0 


11.4 


9.6 


6.6 


1.0 


1.0 


6.7 


0.95 


3.4 


2a 


72.1 


0.7 


7.6 


5.3 


5.5 


2.9 


0.3 


0.4 


4.4 


0.88 


6.1 


2b 


9.02 


1.7 


6.6 


4.7 


4.6 


2.9 


0.1 


0.5 


3.3 


0.92 


2.9 


2c 


1.93 


3.4 


15.8 


11.4 


11.0 


9.4 


0.3 


0.8 


5.4 


0.95 


3.1 


2d 


0.0729 


5.8 


24.6 


17.8 


17.0 


15.5 


0.1 


0.6 


6.8 


0.95 


3.1 


3a 


42.8 


0.6 


7.6 


5.3 


5.4 


2.7 


0.2 


0.6 


4.4 


0.89 


5.2 


3b 


5.45 


1.6 


7.0 


5.0 


4.9 


2.8 


0.2 


0.7 


3.7 


0.92 


2.9 


3c 


1.12 


3.3 


12.2 


8.9 


8.4 


6.0 


0.2 


0.8 


5.6 


0.95 


2.8 


3d 


0.0484 


5.5 


20.5 


15.0 


14.0 


12.0 


0.7 


1.0 


7.0 


0.95 


2.8 


4a 


26.7 


0.8 


7.3 


5.1 


5.2 


2.7 


0.1 


0.6 


4.2 


0.90 


4.3 


4b 


3.68 


1.9 


7.6 


5.4 


5.3 


2.7 


0.1 


0.7 


4.2 


0.93 


2.7 


4c 


0.783 


3.9 


12.6 


9.2 


8.5 


6.1 


0.4 


0.8 


5.7 


0.94 


4.1 


4d 


0.0326 


6.5 


17.2 


13.0 


11.3 


8.9 


0.5 


0.8 


6.8 


0.95 


4.1 


5a 


16.4 


0.7 


7.0 


4.8 


5.0 


2.5 


0.3 


0.5 


4.1 


0.90 


3.7 


5b 


2.21 


1.8 


7.5 


5.4 


5.3 


2.5 


0.1 


0.4 


4.4 


0.94 


3.0 


5c 


0.508 


3.5 


10.2 


7.6 


6.9 


2.8 


0.3 


0.6 


6.1 


0.95 


3.0 


5d 


0.0198 


5.6 


15.0 


11.2 


9.9 


7.1 


0.4 


0.6 


6.7 


0.96 


2.9 


6a 


9.72 


1.0 


7.3 


5.1 


5.2 


3.2 


0.1 


0.6 


3.8 


0.91 


2.6 


6b 


1.65 


2.2 


8.6 


6.2 


5.9 


3.3 


0.2 


0.5 


4.6 


0.95 


2.6 


6c 


0.383 


4.5 


11.1 


8.4 


7.3 


3.7 


0.2 


0.5 


6.1 


0.96 


5.2 


6d 


0.0120 


7.5 


15.0 


11.8 


9.3 


5.4 


0.7 


0.6 


7.3 


0.94 


2.8 


7a 


3.37 


0.8 


6.5 


4.5 


4.7 


2.8 


0.1 


0.7 


3.4 


0.91 


1.5 


7b 


0.614 


1.7 


9.5 


6.7 


6.7 


4.2 


0.2 


0.6 


5.0 


0.96 


1.7 


7c 


0.146 


3.3 


10.6 


7.8 


7.2 


2.8 


0.4 


0.6 


6.4 


0.94 


2.7 


7d 


0.00677 


5.3 


12.0 


9.2 


7.7 


2.5 


0.6 


0.6 


7.1 


0.96 


2.6 



Table 5: Double differential inclusive jet cross sections as a function of Q 2 and P T together with 
their relative errors and hadronisation correction factors. Other details are given in the caption 
to table [3] The bin nomenclature is defined in table [2] 
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2-Jet Cross Section 











total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


cross 


statistical 


total 


uncorrelated 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




section 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




[pb/GcV 3 ] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 a 


10.1 


2.0 


6.7 


4.8 


4.6 


1.9 


0.6 


0.5 


3.9 


0.91 


3.3 


1 b 


5.06 


2.6 


7.6 


5.6 


5.1 


2.2 


0.1 


0.2 


4.4 


0.91 


2.8 


1 c 


1.08 


5.5 


11.5 


8.9 


7.2 


4.6 


0.4 


0.2 


5.4 


0.95 


2.8 


1 d 


0.0435 


8.3 


16.7 


13.1 


10.3 


6.3 


1.0 


0.9 


7.9 


0.96 


2.7 


2a 


6.55 


1.8 


6.7 


4.8 


4.7 


1.7 


0.3 


0.3 


4.1 


0.91 


3.4 


2b 


3.56 


2.4 


8.2 


6.0 


5.7 


3.1 


0.1 


0.6 


4.5 


0.93 


3.0 


2c 


0.767 


5.0 


18.8 


13.7 


12.9 


11.4 


0.2 


0.8 


5.6 


0.95 


2.4 


2d 


0.0289 


7.8 


27.5 


20.2 


18.7 


17.1 


0.2 


0.8 


7.3 


0.94 


3.0 


3a 


3.78 


1.8 


7.4 


5.2 


5.1 


2.1 


0.2 


0.7 


4.4 


0.92 


1.7 


3b 


2.28 


2.2 


7.9 


5.7 


5.4 


2.3 


0.2 


0.7 


4.6 


0.94 


2.8 


3c 


0.441 


4.9 


15.1 


11.2 


10.2 


7.9 


0.3 


0.7 


6.1 


0.94 


2.8 


3d 


0.0201 


7.3 


22.7 


16.8 


15.2 


13.1 


0.7 


1.0 


7.5 


0.96 


3.0 


4a 


2.45 


2.3 


7.3 


5.3 


5.0 


1.8 


0.2 


0.7 


4.4 


0.94 


1.5 


4b 


1.52 


2.7 


8.1 


5.9 


5.5 


2.0 


0.1 


0.8 


4.8 


0.94 


2.2 


4c 


0.327 


5.7 


13.5 


10.3 


8.8 


6.3 


0.4 


0.9 


5.8 


0.92 


1.8 


4d 


0.0135 


8.7 


20.6 


15.8 


13.3 


10.8 


0.6 


0.8 


7.5 


0.97 


3.0 


5a 


1.64 


2.0 


6.8 


4.9 


4.7 


1.8 


0.4 


0.4 


4.1 


0.93 


2.0 


5b 


0.927 


2.5 


8.1 


5.9 


5.6 


2.2 


0.1 


0.5 


4.9 


0.95 


2.2 


5c 


0.213 


5.1 


10.5 


8.2 


6.6 


2.5 


0.3 


0.6 


5.8 


0.96 


2.2 


5d 


0.00763 


7.8 


20.1 


15.2 


13.1 


10.3 


0.5 


0.6 


8.0 


0.96 


2.8 


6a 


0.968 


2.7 


7.2 


5.3 


4.8 


2.3 


0.1 


0.6 


4.0 


0.93 


1.8 


6b 


0.753 


3.1 


8.8 


6.5 


5.9 


2.6 


0.2 


0.5 


5.0 


0.95 


1.9 


6c 


0.138 


6.3 


11.5 


9.2 


6.9 


2.8 


0.4 


0.5 


6.0 


0.95 


4.5 


6d 


0.00521 


9.9 


18.1 


14.5 


10.8 


7.0 


0.7 


0.7 


8.0 


0.97 


4.4 


7a 


0.379 


2.1 


6.7 


4.9 


4.7 


1.7 


0.1 


0.6 


4.0 


0.94 


1.8 


7b 


0.266 


2.4 


8.2 


5.9 


5.6 


1.8 


0.4 


0.6 


5.1 


0.96 


1.8 


7c 


0.0648 


4.8 


11.2 


8.6 


7.3 


1.8 


0.6 


0.6 


6.8 


0.96 


2.7 


7d 


0.00265 


7.0 


15.7 


12.1 


10.0 


5.9 


0.5 


0.5 


7.9 


0.95 


1.8 



Table 6: Double differential 2-jet cross sections as a function of Q 2 and (P T ) of two leading jets 
together with their relative errors and hadronisation correction factors. Other details are given 
in the caption to table [3] The bin nomenclature is defined in table [2] 
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2-Jet Cross Section 















total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 




cross 




statistical 


total 


uncorrelated 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 






section 




uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 






fpb/GeV 2 l 

If / v J 




[%] 


[%1 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 A 




608 




9.7 


35.3 


25.9 


24.1 


23.6 


1.2 


0.6 


4.3 


0.94 


10.2 


1 B 




3450 




3.4 


10.0 


7.4 


6.7 


5.5 


0.5 


0.3 


3.5 


0.96 


7.0 


1 C 




3060 




2.1 


7.2 


5.2 


5.0 


2.8 


0.8 


1.0 


3.6 


0.92 


3.3 


1 D 




687 




3.4 


8.4 


6.4 


5.4 


2.7 


0.4 


1.4 


4.2 


0.87 


3.3 


1 E 




76.6 




7.2 


15.1 


11.8 


9.4 


7.4 


0.6 


1.1 


5.4 


0.85 


3.4 


1 F 




4.40 




15.8 


35.0 


27.2 


22.1 


20.7 


1.3 


0.9 


7.2 


0.92 


3.1 


2 A 




723 




9.7 


12.4 


11.1 


5.6 


4.4 


1.8 


0.8 


2.3 


0.97 


7.6 


2 B 




2470 




3.3 


6.9 


5.3 


4.4 


2.2 


0.7 


0.6 


3.4 


0.97 


7.0 


2 C 




2060 




1.8 


6.8 


4.9 


4.7 


2.2 


0.3 


0.4 


3.9 


0.93 


3.6 


2 D 




441 




3.0 


10.0 


7.3 


6.8 


5.1 


0.3 


0.5 


4.2 


0.89 


1.7 


2 E 




50.7 




6.3 


19.4 


14.4 


13.0 


11.5 


0.9 


0.5 


5.8 


0.87 


1.7 


2 F 




2.31 




15.2 


42.2 


31.7 


27.9 


26.9 


1.7 


0.4 


7.2 


0.87 


2.0 


3 A 




363 




10.4 


13.6 


12.1 


6.3 


1.5 


0.6 


1.0 


5.8 


1.02 


5.2 


3 B 




1430 




3.3 


6.9 


5.3 


4.4 


1.6 


0.9 


0.7 


3.7 


0.97 


4.5 


3 C 




1250 




1.8 


6.7 


4.8 


4.7 


1.7 


0.2 


0.7 


4.1 


0.93 


1.9 


3 D 




261 




3.0 


7.4 


5.5 


4.9 


1.2 


0.2 


0.7 


4.4 


0.90 


2.1 


3 E 




38.2 




5.6 


18.6 


13.7 


12.6 


10.8 


1.4 


0.9 


6.1 


0.88 


2.1 


3 F 




1.21 




16.0 


38.0 


29.2 


24.3 


23.0 


2.3 


1.2 


7.2 


0.87 


4.0 


4 A 




231 




13.9 


15.0 


14.4 


4.2 


1.8 


0.8 


0.6 


3.4 


1.09 


3.0 


4 B 




1020 




4.0 


7.9 


6.2 


5.0 


1.6 


1.0 


0.6 


4.3 


0.97 


5.1 


4 C 




782 




2.2 


6.7 


4.9 


4.6 


1.0 


0.3 


0.7 


4.1 


0.95 


1.9 


4 D 




189 




3.5 


7.7 


5.9 


5.0 


1.1 


0.3 


0.8 


4.5 


0.90 


2.0 


4 E 




26.3 




6.7 


17.2 


13.0 


11.2 


9.2 


1.3 


0.8 


6.0 


0.87 


1.5 


4 F 




0.946 




17.4 


26.7 


22.5 


14.3 


13.0 


1.5 


0.7 


5.5 


0.90 


1.4 


5 A 




82 




14.3 


15.4 


14.8 


4.3 


2.2 


0.6 


0.5 


3.3 


1.04 


4.3 


5 B 




602 




3.7 


7.0 


5.5 


4.3 


2.1 


1.1 


0.4 


3.2 


0.98 


3.8 


5 C 




523 




1.9 


6.8 


4.9 


4.7 


1.6 


0.4 


0.5 


4.1 


0.95 


2.4 


5 D 




120 




3.1 


7.9 


5.9 


5.2 


1.7 


0.4 


0.6 


4.6 


0.90 


1.0 


5E 




16.6 




5.9 


11.2 


8.9 


6.9 


2.9 


1.0 


0.6 


5.9 


0.87 


1.1 


5 F 




0.700 




14.3 


29.5 


23.1 


18.3 


16.4 


1.9 


0.9 


7.7 


0.91 


1.7 


6 A 




33.9 




20.9 


21.6 


21.2 


4.3 


2.1 


1.7 


0.2 


2.9 


0.98 


6.8 


6 B 




431 




5.0 


7.3 


6.2 


3.9 


0.9 


0.7 


0.6 


3.4 


0.99 


2.3 


6C 




335 




2.5 


7.1 


5.3 


4.8 


2.2 


0.4 


0.5 


4.0 


0.95 


2.9 


6 D 




93.8 




4.0 


8.7 


6.7 


5.6 


2.2 


0.3 


0.5 


4.8 


0.91 


0.9 


6E 




9.52 




7.9 


12.6 


10.5 


7.1 


3.0 


1.5 


0.7 


6.0 


0.89 


1.0 


6 F 




0.427 




18.9 


27.7 


23.7 


14.4 


12.3 


1.8 


0.5 


6.9 


0.91 


2.2 


7 A 




4.83 




44.7 


45.5 


45.1 


6.3 


4.8 


0.8 


0.8 


3.5 


1.17 


5.2 


7B 




122 




5.2 


9.6 


7.7 


5.7 


3.7 


0.9 


0.3 


4.0 


0.98 


1.9 


7C 




131 




2.0 


7.2 


5.1 


5.0 


2.9 


0.1 


0.4 


3.8 


0.96 


1.4 


7 D 




36.5 




2.7 


8.8 


6.5 


5.9 


3.1 


0.6 


0.5 


4.8 


0.94 


1.4 


7E 




5.81 




4.9 


11.6 


8.9 


7.5 


3.1 


1.5 


0.4 


6.4 


0.90 


1.5 


7 F 




0.237 




11.8 


22.1 


17.7 


13.2 


10.3 


2.3 


0.4 


7.8 


0.89 


2.7 



Table 7: Double differential 2-jet cross sections as a function of Q 2 and £ together with their 
relative errors and hadronisation correction factors. Other details are given in the caption to 
table [3] The bin nomenclature is defined in tabled 
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3-Jet Cross Section 











total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


cross 


statistical 


total 


uncorrelated 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




section 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 




[pb/GcV 3 ] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 a 


0.780 


4.3 


10.0 


7.6 


6.5 


3.2 


0.3 


0.2 


5.4 


0.78 


11.4 


1 b 


0.951 


3.8 


9.6 


7.2 


6.4 


3.4 


0.3 


0.4 


5.1 


0.81 


6.0 


1 c 


0.240 


7.2 


20.4 


15.2 


13.5 


11.8 


0.7 


0.6 


6.3 


0.82 


5.8 


1 d 


0.0102 


10.9 


24.3 


18.8 


15.4 


12.6 


0.8 


1.1 


8.5 


0.88 


5.5 


II a 


0.321 


4.5 


10.5 


8.0 


6.8 


2.9 


0.6 


0.9 


5.8 


0.78 


8.4 


II b 


0.409 


3.9 


10.0 


7.5 


6.6 


3.1 


0.6 


0.8 


5.5 


0.82 


6.4 


II c 


0.113 


7.3 


17.3 


13.3 


11.2 


8.9 


0.8 


0.9 


6.5 


0.84 


4.8 


lid 


0.00526 


10.3 


24.4 


18.7 


15.7 


13.1 


1.1 


0.9 


8.4 


0.87 


4.7 


III a 


0.126 


5.0 


10.8 


8.3 


6.8 


3.7 


0.5 


0.7 


5.5 


0.76 


10.3 


III b 


0.203 


4.2 


10.9 


8.2 


7.2 


3.9 


0.7 


0.6 


5.8 


0.81 


6.6 


III c 


0.0503 


7.6 


13.9 


11.1 


8.3 


4.0 


0.8 


0.5 


7.0 


0.86 


6.4 


Mid 


0.00231 


10.6 


22.4 


17.5 


14.0 


11.0 


1.1 


0.7 


8.3 


0.87 


4.4 


IV a 


0.0350 


6.8 


13.1 


10.4 


8.0 


5.0 


0.7 


0.5 


5.9 


0.78 


7.1 


IV b 


0.0513 


5.5 


13.9 


10.5 


9.1 


6.6 


1.2 


0.7 


5.9 


0.82 


7.1 


IV c 


0.0167 


9.7 


14.9 


12.5 


8.1 


3.0 


1.3 


0.8 


7.2 


0.85 


5.3 


IV d 


0.000750 


14.1 


20.6 


17.6 


10.6 


7.0 


1.4 


0.7 


7.7 


0.85 


5.3 



Table 8: Double differential 3-jet cross sections as a function of Q 2 and (P?) of two leading jets 
together with their relative errors and hadronisation correction factors. Other details are given 
in the caption to table [3] The bin nomenclature is defined in table [2] 
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3-Jet to 2-Jet Cross Sections Ratio 











total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


3-jet/2-jet 


statistical 


total 


uncorrelated 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




ratio 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 






[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 


0.156 


2.4 


6.8 


5.0 


4.6 


4.2 


0.6 


0.3 


1.2 


0.88 


4.7 


II 


0.162 


2.5 


6.1 


4.6 


4.1 


3.6 


0.7 


0.1 


1.2 


0.88 


5.1 


III 


0.170 


2.4 


5.1 


3.8 


3.4 


2.6 


0.8 


0.1 


1.4 


0.85 


6.9 


IV 


0.151 


3.7 


8.0 


6.2 


5.1 


4.6 


0.8 


0.1 


1.5 


0.85 


6.8 








3-Jet to 2-Jet Cross Sections Ratio 

d(P T ) ' 


dcr 2-jct 
d(P T > 














total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


3-jet/2-jet 


statistical 


total 


uncorrelated 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




ratio 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 






[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


a 


0.0982 


2.4 


4.1 


3.1 


2.6 


1.2 


0.7 


0.1 


1.5 


0.84 


8.5 


b 


0.222 


1.9 


3.2 


2.4 


2.1 


1.2 


0.4 


0.1 


0.7 


0.86 


5.2 


c 


0.276 


3.5 


4.8 


4.1 


2.5 


1.9 


0.4 


0.1 


0.7 


0.89 


3.0 


d 


0.310 


5.1 


6.8 


5.9 


3.3 


2.9 


0.4 


0.1 


0.7 


0.91 


2.9 



3-Jet to 2-Jet Cross Sections Ratio ■fe^y / 











total 


total 


single contributions to correlated uncertainty 


hadronisation 


hadronisation 


bin 


3-jet/2-jet 


statistical 


total 


uncorrelated 


correlated 


model 


electron 


electron 


hadronic 


correction 


correction 




ratio 


uncert. 


uncert. 


uncertainty 


uncert. 


uncert. 


energy scale 


polar angle 


energy scale 


factor 


uncertainty 






[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


1 a 


0.0980 


4.1 


6.7 


5.5 


3.9 


3.2 


0.7 


0.2 


1.4 


0.85 


10.2 


1 b 


0.228 


3.5 


6.6 


5.2 


4.1 


3.7 


0.3 


0.2 


0.7 


0.87 


5.0 


1 c 


0.267 


6.3 


12.0 


9.5 


7.4 


7.2 


0.3 


0.1 


0.8 


0.87 


5.1 


1 d 


0.294 


10.2 


12.4 


11.3 


5.2 


4.8 


0.2 


0.3 


0.9 


0.93 


5.3 


II a 


0.103 


4.4 


6.6 


5.5 


3.6 


2.8 


0.8 


0.2 


1.4 


0.84 


9.2 


II b 


0.215 


3.7 


5.7 


4.7 


3.3 


2.8 


0.4 


0.1 


0.8 


0.88 


6.0 


II c 


0.293 


6.9 


9.7 


8.4 


5.0 


4.7 


0.4 


0.1 


0.5 


0.90 


4.3 


lid 


0.312 


9.2 


12.0 


10.7 


5.5 


5.2 


0.5 


0.2 


0.9 


0.91 


7.6 


Ilia 


0.0966 


4.6 


7.2 


6.0 


4.1 


3.4 


0.9 


0.2 


1.4 


0.81 


11.0 


III b 


0.243 


3.2 


4.9 


4.0 


2.8 


2.2 


0.6 


0.1 


0.9 


0.85 


6.9 


III c 


0.287 


6.9 


9.8 


8.4 


5.1 


4.7 


0.4 


0.1 


1.1 


0.90 


8.4 


III d 


0.359 


9.4 


14.0 


11.9 


7.4 


7.2 


0.5 


0.1 


0.4 


0.90 


8.5 


IV a 


0.0922 


6.7 


10.9 


9.0 


6.2 


5.6 


0.6 


0.2 


2.0 


0.83 


8.8 


IV b 


0.193 


5.5 


11.0 


8.6 


6.8 


6.5 


0.8 


0.1 


0.8 


0.85 


8.8 


IV c 


0.257 


9.0 


13.3 


11.3 


6.9 


6.7 


0.7 


0.2 


0.4 


0.89 


6.6 


IV d 


0.282 


14.1 


16.1 


15.1 


5.6 


5.3 


0.9 


0.2 


0.2 


0.90 


6.3 



Table 9: Ratio of 3-jet to 2-jet cross sections as a function of Q 2 and (P T ) together with their 
relative errors and hadronisation correction factors. Other details are given in the caption to 
table [3j The bin nomenclature is defined in tabled 
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Determination of a s from jets 



Measurement 


a 5 (M z ) 


Uncertainty 


x7 ndf 


experimental 


theory 


PDF 




0.1180 


0.0018 


+0.0122 
-0.0090 


0.0022 


17.5/21 




0.1155 


0.0018 


+0.0121 
-0.0090 


0.0025 


14.3/23 




0.1170 


0.0017 


+0.0090 
-0.0072 


0.0014 


11.0/15 


^jetj C 2-jet ? 0"3-jet 


0.1160 


0.0014 


+0.0093 
-0.0077 


0.0016 


50.6/61 


C3-jet/c2-jet 


0.1215 


0.0032 


+0.0066 
-0.0058 


0.0013 


11.9/13 



Table 10: Values of a s (M z ) obtained from fits to the individual inclusive jet, 2-jet and 3-jet 
double differential cross sections and from a simultaneous fit to all of them and to the ratio 
of 3-jet to 2-jet cross sections. Fitted values are given with experimental, theoretical and PDF 
errors, the normalised % 2 /ndf of the fit is also shown. 
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Inclusive Jet, 2-Jet and 3-Jet Cross Sections 
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Figure 2: Inclusive jet cross sections da^tjdQ 1 
tions da 2 -jet/dQ 2 (c) and da 2 -j e t/d (P T ) (d) and 3-jet cross sections da 3 _j et /dQ 2 (e) and 
dos-jet/d (P T ) (f), compared with NLO QCD predictions corrected for hadronisation. The 
error bars show the total experimental uncertainty, defined as the quadratic sum of the statistical 
and systematic uncertainties. The points are shown at the average values of Q 2 , P T or (P T ) 
within each bin. The NLO QCD predictions are shown together with the theory uncertainties 
associated with the scale uncertainties and the hadronisation (grey band). 
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Inclusive Jet Cross Section 
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Figure 3: Double differential inclusive jet cross sections as a function of Q 2 and Pt, compared 
with NLO QCD predictions corrected for hadronisation. Other details are given in the caption 
to Fig. H 
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2-Jet Cross Section 
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Figure 4: Double differential 2-jet cross sections as a function of Q 2 and (Pt), compared with 
NLO QCD predictions corrected for hadronisation. Other details are given in the caption to 
Fig. El 
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2-Jet Cross Section 
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Figure 5: Double differential 2-jet cross section, as a function of Q 2 and £, compared with NLO 
QCD predictions corrected for hadronisation. Other details are given in the caption to Fig. [2} 
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3-Jet Cross Section 
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Figure 6: Double differential 3-jet cross sections as a function of Q 2 and (Pt), compared with 
NLO QCD predictions corrected for hadronisation. Other details are given in the caption to 
Fig. El 



3-Jet to 2-Jet Ratio 



CM 

P 0.4 



0.3 
0.2 
0.1 
0.0 



H1 


• H1 data 


.Si. 

CNI 


0.4 


H1 


• H1 data 




NLO® hadr 


b 




NLO® hadr 






m 
D 


0.3 




.p 1 T 


* 


-s- -r- 




0.2 












0.1 
0.0 


- • 




(a) 

.... i 






(b) 

. . i ..... 


10 


, ,50 
Q 2 [ GeV 2 ] 






10 <P T )[GeV] 50 



Figure 7: Ratios of 3-jet to 2-jet cross sections as a function of Q 2 (a) and (P T ) integrated over 
the full Q 2 range (b) compared with NLO QCD predictions corrected for hadronisation. Other 
details are given in the caption to Fig. [2l 
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3-Jet to 2-Jet Ratio 
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Figure 8: Ratios of 3-jet to 2-jet cross sections as a function of (Pt) in four different Q 2 ranges 
compared with NLO QCD predictions corrected for hadronisation. Other details are given in 
the caption to Fig. [2} 
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a s from Inclusive Jet Cross Sections 
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Figure 9: Values of a s (Mz) determined using the inclusive jet cross sections measured in 22 
bins in Q 2 and Pt with the A;-factor below 2.5. The error bars denote the total experimental 
uncertainty of each data point. The solid line shows the two loop solution of the renormal- 
isation group equation, a s (Mz), obtained from a simultaneous fit of all 22 measurements of 
the inclusive jet cross sections. The inner band denotes the experimental uncertainty and the 
outer band denotes the squared sum of the PDF uncertainty and the theoretical uncertainties 
associated with the renormalisation and factorisation scales and the model dependence of the 
hadronisation corrections. 
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Figure 10: Values of a s (fi r = a/ (Q 2 + obsV^) extracted from inclusive jet cross section 
(a), 2-jet cross section (b) and 3-jet cross section (c). In each case, the solid lines show the two 
loop solution of the renormalisation group equation obtained by evolving the corresponding 
fitted value of a s (M z ), as summarized in tableQj)! data rows 1-3. Other details are given in the 
caption to Fig.|H 
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Figure 11: Values of a s (fj, r = y (Q 2 + Pj> bJ/2) obtained by a simultaneous fit of all jet 
cross sections in each Q 2 bin (a) and of the ratio of 3-jet cross section to 2-jet cross section (b). 
The solid lines show the two loop solution of the renormalisation group equation obtained by 
evolving the a s obtained from these measurements. For (a) the value of a s is extracted from 
a simultaneous fit of 62 measurements of inclusive jet, 2-jet and 3-jet double differential cross 
sections in bins of Q 2 and Pt ({Pt} f° r 2-jets and 3-jets) with A;-factor below 2.5, see table [101 
4th data row. For (b), a s is extracted from 14 measurements of the 3-jet cross section normalised 
to 2-jet cross section, using only data points with A;-factor below 2.5 for both the 3-jet and 2-jet 
cross sections, see table [T0l 5th data row. Other details are given in the caption to Fig. |9l 
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a s from Jet Cross Sections in DIS 
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Figure 12: Values of a s (fi r = a/ (Q 2 + obs)/2) obtained by a simultaneous fit of all jet cross 
sections in each Q 2 bin of this analysis (squares) together with the fit in different bins at high Q 2 
(circles) [4]. The solid line shows the two loop solution of the renormalisation group equation 
obtained by evolving the a s (Mz) extracted from jets at high Q 2 . Other details are given in the 
caption to Fig.|9] 
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